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Abstract 

We analytically determine the optimal microwave field allowing to reverse the magnetic moment of a 
nanomagnet modeled as a macrospin. This is done by minimizing the total injected energy. The optimal 
field is found to be modulated both in amplitude and frequency, and orthogonal to the magnetic moment 
at any time. Its peculiar shape is directly related to the energy landscape (anisotropy, applied field) and to 
the damping parameter. For typical values of the damping parameter, a weak microwave field is sufficient 
to induce switching. The results are in good agreement with the fields previously obtained using the 
optimal control theory. The total injected energy is found to be proportional to the static energy barrier 
between the initial state and the saddle point and to the damping parameter. This result may be used as 
a means for probing the damping parameter and surface effects in real nanoparticles. 

PACS numbers: 75.10.Hk, 75.75.Jn, 84.40.-x 
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I. INTRODUCTION 



Magnetic recording is a key technology in the field of high density information storage. In order 
to increase thermal stability, small nanoparticles with high anisotropy may be used. However, high 
fields are then needed to reverse the magnetization but these are difficult to achieve in current 
devices. To overcome this so called magnetic recording trilemma, several solutions are being 
proposed. The most investigated route, and which already leads to industrial applications, is the 
heat-assisted magnetic recording 1 . It consists in heating the particles by a laser which decreases 
the energy barrier between the two energy minima and thereby the switching fields. However, to 
avoid a loss of information, the heating must be very localized and followed by a very fast cooling, 
and as such these devices must be coupled to powerful heat dissipation systems. 

An alternative solution is to assist the switching by a microwave (MW) field. In 2003 Thirion et 
al.— showed that the combination of a DC applied field (static field) well below the switching field 
with a small MW field pulse can reverse the magnetization of a nanoparticle. Further experimental 
and theoretical studies-- proved that this technique is more efficient if the frequency of the MW 
field matches the ferromagnetic resonance frequency of the nanoelements. The issue of state 
transition assisted by a MW field has previously been addressed in many areas of physics and 
chemistry, especially in the context of atomic physics. Furthermore, Meerson and Friedland 6 
studied the ionization of Rydberg atoms and predicted that a quasi-monochromatic oscillating 
electric field, whose frequency slowly decreases with time, should be far more efficient than a purely 
monochromatic field for achieving ionization. Measurements^ and both analytical and numerical 
calculations^ have confirmed this fact for nanomagnets. 

In our previous work-, we described a numerical method based on optimal control theory which 
renders an exact solution for the MW field that is necessary for the switching of a nanomag- 
net within a given potential energy. The formulation of this method consists in defining a cost 
functional and minimizing it using the conjugate gradient technique. Our results confirmed that 
a weak MW field, modulated both in amplitude and frequency, can induce the switching of the 
magnetization. Furthermore, the injected energy has been found to increase with damping. 

The aim of the present study is to provide analytical expressions and to compare them with our 
numerical results, by using simple energy considerations. Moreover, the analytical developments 
presented here confirm the effects observed numerically and provide clear interpretations for the 
underlying physical processes. In the first part, we analytically determine the optimal MW field 
and show that it is directly related to the energy landscape (anisotropy, applied field) and to the 
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damping parameter. We then investigate the trajectory of the magnetization in the presence of 
the optimal field and show that it can be described by the Landau-Lifshitz-Gilbert equation with 
a negative damping parameter. In the second part, the analytical results are compared directly 
with the results obtained numerically using the optimal control theory. 

II. ANALYTICAL CALCULATION OF THE OPTIMAL MICROWAVE FIELD 

We consider a nanomagnet with homogeneous magnetization which can be modeled by a vector 
M = M s m, where M s is the saturation magnetization and ||m|| = 1. This nanomagnet is 
characterized by a given anisotropy (uniaxial, biaxial, cubic.) and a damping parameter a. In 
the presence of a static magnetic field H lower than the Stoner-Wohlfarth switching field, the 
potential energy surface presents several minima separated by saddle points. 

At the initial time ti, we assume that the magnetization is in a minimum Mj = M^nij. Adding a 
microwave (MW) field B(t) can then induce switching to another (final) minimum My = M s u\f. 
Our aim is to find the optimal field B°P^(t) that achieves switching in a given time tf while 
minimizing the energy injected into the system. 

For the sake of simplicity, we introduce the reduced fields h = H/iJan and h(t) = B(t)//Jan, 
where i^an = 2K/fi M s is the anisotropy field and K the anisotropy constant of the nanomagnet. 
We will also define the reduced time r = 7-f/an^, where 7 = 1.76 x 10 11 (T.s) -1 is the gyromagnetic 
factor. For instance, for a cobalt particle of 3 nm diameter with K 2.2 x 10 5 J.m -3 and 
M s ^ 1.44 x 10 6 A.m' 1 we have n H m « 305 mT and t/r w 1.86 x KT 11 s. 

A. Energy and time trajectory in the presence of a microwave field 

If only the static field is applied, the energy density of the system reads £o(m, h) = £an( m ) — 
m • h where £an is the anisotropy energy density. The reduced effective field is then defined by 
h e g = — d£ /dm. If we add a MW field, the energy density becomes 



£(m, h, b) = £an(m) - m • (h + b(r)) 



(1) 



and the reduced total effective field now reads 



C(r) 



dm 



h eff + b(r). 



(2) 
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The time trajectory of the magnetization can be described by the driven Landau-Lifshitz-Gilbert 
equation 

I 1 + tt 2 ) —j— = -m x C(r) - am x (m x £(r)) . (3) 

This allows us to express the energy variation of the system as follows 

dS . , s dm db , dm c? , , s , . 

Next, we define the mobile frame (m, u, v) attached to the magnetization with u = T/T and 
v = m x T/T, where T = m x h g g and T = ||m x h e g||. The MW field can then be decomposed 
as b(r) = 6 m (r)m + 6 u (r)u + b v {r)\. In this frame, Eqs. ([3]) and (jlj) respectively become 

, 2 ,dm ( ab u + b v \ ( b u -ab v \ 

( 1 + «)^7=(- 1 + — f — j 1" — t — H mxT )' ( 5 ) 



d£_ -aT -b u + ab v db m 

dr ~ 1 + a 2 dr' U 

We note that the parallel component of the MW field b m has no direct effect on the magnetization 
trajectory and that its time derivative appears only in the energy variation. 



B. Optimization of the MW field 

In order to find the optimal MW field fulfilling the requirements described earlier we proceed in 
two steps. First, we define the critical MW field which allows us to maintain the precession of the 
magnetization by compensating the effects of damping (sec. Ill B ip . This field represents the lower 
limit for the optimal field sought. Using this result, we find the optimal MW field minimizing 
the injected energy (sec. Ill B 21) and check that it can induce switching of the magnetization (sec. 
HTB3]) . 



1. MW field maintaining the precession 

In order to induce switching the MW field must at least compensate for the effect of damping, 
which tends to take the magnetization back to the initial equilibrium position. If the compensation 
is complete the energy variation of the system dS/dt vanishes at any time thus reflecting the 
conservation of energy. According to Eq. ([6]) an infinity of MW fields leads to a full compensation 



4 



of damping. For instance, any field that is orthogonal to the magnetization so that b m = and 
satisfying the equation — b u + ab v = —aT will do. 

Among these MW fields the critical field can be defined as the one that minimizes the power 
injected in the system. The latter is proportional to the squared intensity of the MW field, i.e. 
Pi(r) = b 2 (r) = b 2 m (r) + 6„(r) + b 2 (r). Using the method of Lagrange multipliers we define the 
functional 

L [b m (r), b u (r), b v (r), A(r)] = p 2 (r) - A(r)^. (7) 

Assuming that db m (r)/dr does not depend explicitely on b m (r), the minimization of this func- 
tional leads to 



bjr) = 0, 



a 

a 2 



W = -rhi^ (8) 



The critical MW field thus reads 



b crit /x = [_ T + a m x T ] ( 9 ) 

1 + a 2 



and the injected power is then 

pf t (r) = -^— 2 T 2 . (10) 
1 + a z 



In the presence of this MW field the magnetization precesses around the equilibrium position 
with a constant angle. This critical field represents a lower limit. Indeed, if the injected power is 
smaller than pf ri ^(r) the magnetization goes back to the initial equilibrium position. 



2. MW field minimizing the total injected energy 

In order to minimize the total injected energy we have to make a few preliminary assumptions 
concerning the shape of the MW field. Considering the result of the previous section we limit our 
search to the family of MW fields defined by 

b m {r) = f3 m T 
bu(r) = P U T 
b v (r) = (3 V T 
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where /3 m , (3 U and /3 V are constant parameters. In the presence of such a MW field the energy 
variation reads 



d£ 



-a- (3 U + a(3 v 



rp2 _ 



1 + a 2 ' dr 

The total energy injected to the system can be defined as E = j T } pi(r)dt. Therefore, 



E 



{fa + + ft) T 2 dr 

J Tj 



(12) 



rf 



(fin + + ft) 



l + a 2 



dE db r 



—a — (3 u + a(3 v \dr dr 



dr 



[£(r f ) - Sin) + b m (r f ) - b m (ri)\ ■ 
-a - p u + ap v 

Hence, the energy is minimal if (3 m = 0, (3 U = —2a/ (1 + a 2 ) and (3 V = 2a 2 / (1 + a 2 ). Conse- 
quently, the optimal MW field is 

2a 



[-T + amx T] = 2b crit (r). 



(13) 



l + a 2 

The optimal field is twice the critical field determined previously, see Eq. It is orthogonal 
to the magnetization at any time and its magnitude reads 

2a 



b °Pt( r ) 



VTT^ 2 



T. 



(14) 



The total injected energy is then 



E = 4a [8 (r f ) - £ (r,)] 



(15) 



According to these results, both the optimal field magnitude and the total injected energy 
increase with damping. This confirms the fact that the MW field must compensate for the effects 
of damping so as to induce switching. 



3. Trajectory of the magnetization in presence of the optimal MW field 

In order to check whether the optimal MW field obtained in the previous section induces switch- 
ing of the magnetization as required, we now investigate the time trajectory of the magnetization. 
In the presence of this field, Eqs. (j5]) and ([6]) respectively become 

, 9 x dm , , 

(l + a 2 ) — = -T + am x T, 16 

v ' dr 
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dS 
d^ 



1 + a 



a 




(17) 



The first equation is similar to the Landau-Lifshitz-Gilbert equation but with a negative damp- 
ing parameter: it describes an amplified precession. The precession frequency is equal to the proper 
frequency of the magnetization. At any time the MW field is proportional to the derivative of the 
magnetization: b(t) = 2adm/dt. This is in agreement with the results of Sun et al. 4 . 



T = 0. Therefore, both the derivative of the magnetization and the MW field vanish. This has 
two consequences: (i) small fluctuations around the equilibrium position are necessary to start the 
amplified precession and (ii) the magnetization cannot cross a saddle point. The optimal MW field 
thus induces the motion of the magnetization from an initial state m^, close to an energy minimum, 
to a final state mj close to a saddle point. Fluctuations then allow the system to cross the saddle 
point and afterwards the damping takes up to lead the magnetization down to the second energy 
minimum. If the energy landscape is complex with several barriers successive pulses may then be 
necessary to induce switching. 

At both the initial and the final states the MW field is close to zero. The difference £(t/) —8 fa) 
is thus close to the static energy barrier between the saddle point and the initial state ASq = 
8 (rf) — £o( r i) an d Eq. (fT5|) becomes 



The total injected energy is therefore proportional to the energy barrier to be overcome. Hence, 
if the static field is close to the switching field, a very weak MW field can induce switching. 

C. Uniaxial anisotropy and longitudinal static field 

In this section we study the trajectory of the magnetization in the presence of the optimal MW 
field for a nanoparticle with uniaxial anisotropy and a longitudinal static field. 

We consider a nanomagnet with uniaxial anisotropy with easy axis in the z direction. The 
anisotropy energy density is then £an( m z) — —km 2 z /2 with k — 1. The static field is applied 
in the (—z) direction with a magnitude < h < k. The magnetization is initially close to the 
metastable minimum and its z component is thus m = m z {ri = 0) ~ 1. Setting x = h/k, the 
static energy of the system is 



At the minima and saddle points the effective field h e g is parallel to the magnetization so that 



E = 4a AS, 





(19) 
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The saddle point then corresponds to m z = x, so the static energy barrier between the latter 
and the initial metastable state is A£ = k (1 — x) 2 /2. 

For this system, the effective field reads h g g = k {—x + m z ) e z . Projecting Eq. f fl~6l) onto the z 
axis then yields 



dm z 



a 



m x TVe, 



a k 



dr 1 + a 2 " ' 1 + a 2 

In order to simplify the expressions we introduce the integral 

-du 1 



(— x + m z ) (l - mf) . 



(20) 



l\m> 



In 



1 + m z f x (1 - m z 



(—x + m z y 



(~x + u)(l-u 2 ) 2(l-x 2 ) 

Solving Eq. ([20]) with the initial condition m z (t = 0) = m leads to the equation 

ak 



+ C 



te 



(21) 



Km, 



I (ma) 



T. 



(22) 



1 + a 2 

This equation can be analytically solved for m z only if x = (no static field). Otherwise, the 
evolution of m z with time can be obtained numerically, see Fig. [H As predicted previously, for 
long times the magnetization goes towards to the saddle point but never reaches it since I(m z ) 
diverges for m z = x. 

>From Eq. (I14p we can express the optimal MW field intensity in terms of m z as follows 

2ak 



m, 



VT+a- 



-.(-x + m z )y / l 



ml- 



(23) 



The time evolution of the MW field intensity is plotted in Fig. [TJ We note that the pulses 
follow neither a Gaussian nor a Lorentzian function. The peak intensity is reached for m z = 
| (x + y/8 + x 2 ) and is given by 



■ opt 

J max 



ak 



( , \ / (x + y/8 + X 2 ) 

{-3x + V8 + ^j W 1 - ^ ^— >-. 



(24) 



2^/l + a 2 

>From this, we can see that the peak intensity &nrax decreases with x (see Fig. |6]). Indeed, for 
higher magnitudes of the static field, the energy barrier A£ between the metastable state and the 
saddle point is lower, so that a lower energy is needed to reach the saddle point. Since < x < 1 
the peak intensity is limited as follows 

y/2ak 



6 opt 
"max ^~ 



VT+a 2 



(25) 
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Figure 1: Optimal MW field intensity 



b°P (t) (upper panel) and trajectory of the magnetization m z (r) 



(lower panel) for a longitudinal static field with magnitude x = h/k. Parameters: a = 0.05, tuq = 0.99998. 



Hence, for low values of the damping parameter a, the intensity of the optimal MW field is 
small. This fully confirms the results of our numerical study 9 . Using Eqs. ( l2TT) and f l23l) . we can 
also obtain analytically the peak duration At, defined as the full width at half maximum 

1 + a 2 



ak 



-9 (x) , 



(26) 



where g (x) is a cumbersome function of x. This characteristic time increases with x (see Fig. [6]). 
For high values of x, the switching will thus require a very low field but a very long time, as can 
be seen in Fig. [TJ Moreover, the characteristic time decreases with damping. 



The area below the curves 
A = 



t=oo 



4=0 



b opt ( 

b °Pt( 7 



is 



dr 



m z =mo 



dr 

dm. 



dm z 



(27) 



2Vl + a 2 [arccos (x) — arccos (mo)] = 2Vl + a 2 (Of — 0{ 



where 0i and Of are the polar angles of the magnetization respectively at the initial state and 
saddle point. This area is thus proportional to the "angular distance" that the magnetization must 
cross to reach the saddle point. For increasing values of the static field magnitude x, this area 
decreases since the saddle point comes closer to the initial state. 
The z component of b°P^ , given by Eq. (fTBl) . is 

2a 2 k 



1 -1 

>From Eq. (123]) we can see that the ratio 



{—x + m z ) (l — m^) . 



(28) 



■ opt, 



b opt, 



m. 



is 



m 



?), whose 



upper limit is ^j=?- Consequently, for small values of the damping parameter a, the component 
of the time- dependent field along the anisotropy easy axis can be neglected and the optimal field 
lies in the xy plane. 

We now define the precession phase of the magnetization as ip(r) = arctan (m y (r) / m x (T)) . 
Projecting Eq. ( JT6|) on the x and y axes leads to the relation 



dip 
dr 



^ffl 



k 



:(-x + m z ). 



(29) 



1 + a 2 1 + a 2 

This precession frequency is equal to the proper frequency of the magnetization, obtained by 
solving Eq. [3] in the absence of a MW field. At the initial state, the precession frequency is close 
to the FMR frequency ojfmr — k(l — x) / (1 + a 2 ). It then decreases towards zero, following the 
curvature of the energy well. 

For small values of a, the optimal field lies in the xy plane and its phase can be defined by 
tp(t) = arctan (by^ (t)/bx^ (t)j. It can be shown that 



tan<£>(r) 



m x {r) + am y {T) m z (r) 



m x {T 



COt 



(30) 



-m y (r ) + am x (r ) m z (r) m y (r) 
This implies that the time-dependent field and the magnetization are synchronized with (p(r) ~ 
(p(r) + 7r/2. Hence, the frequency of the time-dependent field is equal to the proper precession 
frequency of the magnetization. 
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Finally, using Eqs. (120]) and (123]) . the total injected energy can be computed. As shown 
previously in Eq. (ITS]) , it is proportional to the damping parameter and to the energy barrier A£ . 



E 



r+°° , 2 

/ h ^ (r) dr = 2ak (x - l) 2 = 4aA£ . 
Jo 



(31) 



III. COMPARISON WITH THE NUMERICAL RESULTS 

In order to check the validity of our model, we have carried out numerical calculations based 
on the optimal control theory, as presented irA Our model system is a particle with uniaxial 
anisotropy along the z axis with k = 1. The numerical method consists in minimizing the cost 
functional 

^[m(T) l b(r)] = i||m(r / )-m / || 2 + | J drh\r) 

o 

along the trajectory given by the Landau-Lifshitz-Gilbert equation (j3J), with mj is the target 
magnetization (stable minimum), va{r f) is the magnetization reached at time Tf and rj is a nu- 
merical control parameter. The numerical problem is solved using the modified conjugate gradient 
technique supplemented by a Metropolis algorithm. 

Unless otherwise specified, the numerical parameters used in the current study are: initial 
reduced time n = 0, final reduced time Tf — 800 (corresponding to a few ns in real time), number 
of points iV = 15000, so the sampling time (t/ — r»)/(iV — 1) is about 0.05; damping parameter 
a = 0.05; control parameter rj = 0.01. For a better comparison with the analytical results, the 
MW field was allowed to move in any direction during the optimization. 



A. Reference calculation 



A first numerical optimization was carried out for a static field with magnitude h = 0.5 (x = 0.5) 
applied in the (—z) direction. The results are in good agreement with the analytical calculations 
of part Ull As can be seen in Fig. [2j the optimal MW field is modulated both in amplitude and 
frequency. The MW field is mainly in the xy plane and its magnitude is small (less than 0.03, 
corresponding to a few mT in real field) as expected since the damping parameter is small. The 
pulse starts at about r = 350 and vanishes as soon as the magnetization crosses the saddle point 
(at about r = 670). Next, the damping takes up to lead the magnetization to the more stable 
energy minimum. 
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Fig. [3] shows that the MW field intensity obtained numerically is in good agreement with the 
analytical result in Eq. ([2"3"[) . From r ~ 570, the MW field intensity is slightly higher numerically, 
which induces m z to decrease faster and finally the magnetization to cross the saddle point. Fig. H] 
confirms that the magnetization precession and the MW field are synchronized, the initial frequency 
being close to the FMR frequency. The time evolution of the frequency is similar to the evolution 
of m z (Figj3]), since both values are related by Eq. (j29|) . Consequently, after r f» 570 the numerical 
frequency is lower that the analytical frequency. 



B. Effect of the static field magnitude and direction 



The MW field b(r) has been optimized numerically for several magnitudes and orientations of 

the static field h. For each configuration, the energy barrier A£ Q and the total injected energy 

, 2 

E = p h°P l {r) drhave been computed numerically, and are reported in Fig. [SJ As predicted 
analytically by Eq. ([TBI) , the injected energy is found to be proportional to the energy barrier and 
to 4a. 

In the case of a static field applied along (—z), the shape of the pulse can be directly compared 
with the analytical results of sec. Ill C| see Fig. |6j The numerical and analytical results are in 
good agreement. As expected, the maximal peak intensity fe^fax decreases and the pulse duration 
At increases rapidly when the magnitude of the static field x increases. 



C. Effect of damping 

As predicted by Eqs. fl!8j) and (I31D . for a given static energy barrier ASq, the injected energy 
is proportional to the damping parameter a. This has been checked numerically by varying the 
damping parameter from 0.015 to 0.30 (Fig. [7]). In these calculations, the static field h is applied 
in the (—z) direction with the magnitude x = 0.5. 

Fig. [HJ shows that for low values of a, the pulses height decreases but their duration increases, 
in agreement with Eqs. ([24"j) and (1261) . For an undamped system, the switching should thus be 
infinitely long, so our analytical and numerical methods are not adapted to describe such a system. 

As can be observed, for very low values of the damping parameter a, a discrepancy between the 
analytical calculations and the numerical optimization is observed. Indeed, since the optimal peak 
duration becomes very long, the number of numerical points iV must be increased, so the conjugate 
gradient algorithm becomes less efficient. However, as can be seen in Fig. [7[ this discrepancy has 
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Figure 2: Numerical results for the reference calculation: optimal MW field (upper panel) and magneti- 
zation time trajectory (lower panel). 

a negligeable effect on the injected energy. 
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Figure 3: Comparison between the analytical and numerical results for the reference calculation: optimal 
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(higher panel) and z component of the magnetization (lower panel). 



IV. CONCLUSION 



We determined analytically the optimal microwave field allowing for the switching of the mag- 
netization of a monodomain nanoparticle with uniaxial anisotropy while minimizing the injected 
energy. This study provides a clear interpretation of the results obtained numerically using the 
optimal control theory-, especially the simple dependence of the pulse on the damping parameter. 

Our results show that the optimal MW field is modulated both in amplitude and frequency, 
and is directly proportional to the derivative of the magnetization. For typical values of the 
damping parameter (a < 1), a weak MW field is sufficient to induce switching. Its role is to 
drive the magnetization from the initial state towards the saddle point; the time trajectory of the 
magnetization can then be described as an amplified precession. After crossing the saddle point 
the damping induces the relaxation to the final state. The injected energy is then proportional to 
the damping parameter and to the static energy barrier between the initial state and the saddle 
point. 

For a nanomagnet with uniaxial anisotropy placed in a longitudinal static field, the shape of 
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Figure 4: Comparison between the analytical and numerical results for the reference calculation: precession 
frequencies of the magnetization and MW field. 

the MW field pulse has been obtained analytically. We have shown that the optimal MW field 
pulse becomes smaller but much more spread when the damping decreases. 

In the case of more complex energy landscapes (biaxial or cubic anisotropics), the switching 
is likely to be triggered by a succession of MW field pulses. This hypothesis will be later tested 
numerically. This study could then be extended to small nanomagnets where surface effects can 
not be neglected, using the effective one-spin model (EOSP)— ~— . 

The MW fields that we optimized have an amplitude and a frequency which vary slowly, and 
can be reproduced experimentally using a frequency generator. Consequently, our theoretical 
results could be used to probe the damping parameter and assess the role of surface effects in real 
nanoparticles. The dependence of the MW field on the energy landscape might be used to address 
directly a given nanoparticle in a polydisperse assembly. 
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Figure 7: Total injected energy E with respect to the damping parameter a. 
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Figure 8: Optimal MW field pulse for several values of the damping parameter a. For a = 0.015, the final 
time has been increased (rt = 1600) without changing the sampling time. 
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